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Abstract 

We studied magnetism of a spin-1 substance Li2Ni2Mo30i2. The spin system consists of distorted 
honeycomb lattices and hnear chains of Ni^"*" spins. Li^ ions enter about 25 % and 50 % of 
honeycomb and chain Ni sites, respectively, creating disorder in both the spin subsystems. A 
magnetic phase transition occurs at Tc = 8.0 K in the zero magnetic field. In low magnetic fields, 
the magnetization increases rapidly below T^, decreases below 7 K, and finally becomes negative at 
low temperatures. We determined the magnetic structure using neutron powder diffraction results. 
The honeycomb lattices and linear chains show antiferromagnetic and ferromagnetic long-range 
order, respectively. We investigated static and dynamic magnetic properties using the local probe 
technique muon spin relaxation. We discuss the origin of the negative magnetization. 
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I. INTRODUCTION 



Several antiferromagnets having plural magnetic ions or magnetic-ion sites show an in- 
teresting magnetic long-range order (LRO). In Cu2Fe2Ge40i3, Cu^"*" spin dimers (spin 1/2) 



are coupled to Fe^^ spin chains (spin 5/2). [l|, |2| An indirect Fe-Fe exchange coupling via 



m 

Cu dimers was observed. [3i] This result reveals that the Cu dimers play the role of nonmag- 
netic media in the indirect magnetic interaction. Dynamics of dimer excitation in staggered 
and random fields generated by Fe spins are also experimentally identified. [^] In CU3M02O9, 
three crystallographically independent Cu^"*" sites exist (Cul, Cu2, and Cu3). Spins (1/2) 
on Cul sites form antiferromagnetic (AFM) chains, [sl Two spins on neighboring Cu2 and 
Cu3 sites form an AFM dimer. Each AFM chain is coupled to AFM dimers. Only a com- 
ponent parallel to the b axis of magnetic moments on Cul sites is ordered below the AFM 
transition temperature Tn = 7.9 K in the zero magnetic field, [s] Perpendicular components 
are ordered in magnetic fields parallel to the a or c axis. A canted AFM LRO is stabilized 
only in finite magnetic fields. It is inferred that magnetic competition causes this unique 
LRO. In addition, the canted AFM LRO disappears in CU3M02O9 with a very small amount 
doping {e.g., 0.2 % Zn for Cu sites). 

Negative magnetization is an interesting phenomena. It was first observed in 
Lio.5(FeCr)2.504 spinels [3] and was explained using the Neel model, jsl In this model, two 
collinear sublattice magnetizations (Ma and Mb) with different temperature (T) dependence 
are considered. |Ma| is larger than |Mb| below Tn, while |Mb| is larger than |Ma| below 
a compensation temperature Tcomp- If directions of the sublattice magnetizations are fixed, 
magnetization is observed as a negative value below Tcomp- Several compounds showing 
negative magnetization have been found. Some of them are considered to be explainable by 



the Neel model, 



18l | while for the other compounds the origin is different or has not been 



determined explicitly. 19|-|37| Details of the origins will be described later. 

We have paid our attention to Li2Ni2Mo30i2. Two crystallographically independent Ni^"*" 



(spin 1) sites exist (Ml and M2) as represented in Fig. 1. 2 



The Ml and M2 sites form 



distorted honeycomb lattices and linear chains, respectively. About 25 % of Ml sites and 
about 50 % of M2 sites are occupied by non-magnetic Li"*" ions as shown by the formula 
Li(Lio.5Nii.5)(Lio.5Nio.5)Mo30i2. Unexpectedly, we found negative magnetization at low T. 
To investigate the origin of the negative magnetization, we determined the magnetic struc- 
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FIG. 1: (Color online) A schematic drawing of Ml {8d) and M2 (4c) sites in Li2Ni2Mo30i2. The 
space group is orthorhombic Pnma (No. 62).j38] These sites are occupied by Ni^^ (spin 1) or 
Li^" (spin 0) ions. The occupation ratio of Ni^+ ions is 0.757 at Ml and 0.566 at M2. The ratio 
between Ni numbers at Ml and M2 sites is 0.728 : 0.272. Three short Ml-Ml bonds exist and 
form distorted honeycomb lattices. The Ml-Ml distances are 3.03, 3.08, and 3.75 A at 10 K for the 
first, second, and third shortest Ml-Ml bonds, respectively. The respective exchange interaction 
parameters are defined as Ji, J2, and J3. The shortest M2-M2 bonds form linear chains. The 
M2-M2 distance is 2.53 A. Distances in the other Mi — Mi bonds (i = 1 or 2) are larger than 5.06 
A. Distances in M1-M2 bonds are larger than 5.19 A. 



ture of Li2Ni2Mo30i2 using the neutron powder diffraction technique and determined its 
static and dynamic magnetic properties using the local probe technique muon spin relax- 
ation (/iSR). 



II. EXPERIMENTAL METHODS 

Crystalline powder of Li2Ni2Mo30i2 was synthesized using a solid-state-reaction method 
at 973 K in air for 144 h with intermediate grindings. We used an isotope ^Li (purity: 
99 %) for neutron-diffraction samples to decrease absorption of neutrons. We confirmed 
formation of Li2Ni2Mo30i2 using an X-ray diffractometer (JDX-3500; JEOL). Li2Ni2Mo30i2 
is insulating at room temperature. 

We measured magnetization using a superconducting quantum interference device 
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(SQUID) magnetometer (MPMS-5S; Quantum Design). We entered powder in paraffin 
molten by heating in the zero magnetic field and fixed finally the powder in solid paraffin. 
We measured magnetization of the powder embedded in solid paraffin. Probably, powder 
directions are randomly distributed and powder reorientation does not occur in the magneti- 
zation measurements. Specific heat was measured using relaxation technique with Physical 
Property Measurement System (PPMS; Quantum Design). We used a sintered pellet in the 
specific-heat measurements. 

We determined the magnetic structure of ''Li2Ni2Mo30i2 from neutron powder diffraction 
data. The experiments were conducted using the high-resolution powder diffractometer for 
thermal neutrons HRPT 39 1 (wavelength A = 1.886 A) and the high- intensity cold neutron 



powder diffractometer DMC (A = 4.206 A) at the Swiss spallation neutron source SINQ in 
PSl. Powder was filled into a vanadium container with 8 mm diameter and 55 mm height. 
Rietveld refinements of diffraction data were performed using the FULLPROF Suite program 
package. j4ol| Symmetry analyses of possible magnetic configurations were conducted using 
the program BASIREP in the FULLPROF Suite program package. 

The muon spin relaxation (/iSR) measurements were performed using the Dolly spec- 
trometer at the ttEI beamline at the Swiss Muon Source (SyuS) in PSl. In a yuSR experiment 
nearly 100 % spin-polarized muons are implanted into the sample one at a time. The pos- 
itively charged fi~^ thermalize at interstitial lattice sites, where they act as magnetic micro 
probes. In a magnetic material the muon spin precesses about the local magnetic field 
B at the muon site with the Larmor frequency = '~f^/ (2,71)3 (muon gyromagnetic ratio 
7^/27r = 135.5 MHz T^^). With a lifetime of = 2.2 /is the muon decays into two neutrinos 
and a positron, the latter being predominantly emitted along the direction of the muon spin 
at the moment of the decay. Measurement of both the direction of positron emission as 
well as the time between muon implantation and positron detection for an ensemble of sev- 
eral millions of muons provides the time evolution of the muon spin polarization P{t) along 
the initial muon spin direction. In a powder, 2/3 of the local magnetic field components 
are perpendicular to the yU+ spin and cause a precession, while the 1/3 longitudinal field 
components do not. The damping of the oscillation is a measure of the width of the static 
field distribution experienced by the muon ensemble. In a static magnetic environment, the 
1/3 fraction of the muon spin polarization is conserved. In a dynamic magnetic materials 
however also this so-called 1/3-tail relaxes and the relaxation rate can be, in certain limits. 
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FIG. 2: (Color online) Temperature dependence of the specific heat of Li2Ni2Mo30i2 in the zero 
magnetic field. 



related to magnetic fluctuation rates. Since /xSR is a local probe technique, the amplitudes 
of the different signals observed are a measure of the corresponding volume fractions. With 
this, it provides a direct measure of magnetic volumes which is not easily possible with scat- 
tering techniques. For further details of the /iSR technique, we refer to the recent review in 



Ref. 



4l|. 



III. EXPERIMENTAL RESULTS 



Figure 2 portrays the T dependence of the specific heat in the zero magnetic field. A peak 
is apparent at 8.0 K. As it will be shown later, magnetization increases rapidly on cooling 
below 8.0 K. Therefore, the peak indicates occurrence of a magnetic phase transition at 
Tc = 8.0 K. 

Figure 3 (a) shows the T dependence of the magnetization M of Li2Ni2Mo30i2 in a very 
weak residual magnetic field. The value of the magnetic field was estimated as 6.3 x 10"^ T 
(0.63 Oe) from paramagnetic magnetization values in this field and 0.1 T. The magnetization 
increases rapidly below 8.2 K, shows a maximum around 7.2 K, and decreases below 7.2 K. 
It should be noted that the magnetization is negative below 4.2 K. We did not observe 
any difference between magnetizations measured for increasing and decreasing T. Figure 4 
represents the magnetic field H dependence of M at 1.7 K after the residual field cooling 
process. We can see a small residual magnetization and hysteresis. The rapid increase 
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below 8.2 K and small residual magnetization indicate appearance of a canted AFM LRO 
or ferromagnetic (FM) LRO with small magnetic moments. Figure 3 (b) shows the T 
dependence of M in the field of 1 x 10^^ T. We can see a hysteresis between magnetizations 
measured in the residual field cooling (RFC) and field cooling (FC) processes. The absolute 
value of M in the FC process is larger than that in the RFC process. In 0.1 T [Fig. 3(d)], 
only the FC magnetization becomes negative. In 0.2 T [Fig. 3(e)], both the RFC and FC 
magnetizations are positive. In 1 T [Fig. 3(f)], the hysteresis is small. Figure 3 (g) shows 
the T dependence of paramagnetic susceptibility in 0.1 T. The susceptibility above 20 K 
obeys the Curie- Weiss law. We determined the value of the Curie constant as 1.12 emu 
K/mol Ni from the data above 200 K. The g value was calculated as 2.11 on the assumption 
that the spin value is 1. This g value is reasonable for Ni^+ spins. 

Figure 5 depicts the neutron powder diffraction pattern of paramagnetic ''Li2Ni2Mo30i2 
recorded using the HRPT diffractometer with A = 1.886 A at 10 K, which is slightly higher 
than Tc = 8.0 K. The refinement based on the crystal structure of Li2Ni2Mo30i2 as deter- 
mined by room temperature single-crystal X-ray diffraction [38] well fits the experimental 
neutron diffraction pattern at 10 K. Structural parameters are presented in Table I. 

Figure 6(a) shows difference between two neutron powder diffraction patterns of 
^Li2Ni2Mo30i2 at 2.0 K and 10 K collected using the DMC diffractometer (A = 4.206 
A). The reflections of the difference pattern can be indexed in the chemical cell with the 
propagation vector k = 0. The proflle matching le Bail flt, in which the integral intensities 
are the reflned parameters, shows that all the peaks are well described with k = with the 
reliability factors -R„p = 2.34, i?exp = 1.72, and = 1-86. The new reflections appear below 
Tc. Therefore, they must be magnetic reflections. 

Using the determined propagation vector, we performed a symmetry analysis according 
to Izyumov et al. 42(] to derive possible magnetic conflgurations for the space group Pnma. 
The space group has eight one dimensional irreducible representations (IRs). The observed 
magnetic patterns were compared with calculated patterns using the structural parameters 
determined from the structural reflnement. After sorting out the basis functions of all 
eight IRs, we found that only well flts the experimental pattern. The IR has the 
following characters 1, -1, 1, -1, 1, -1, 1, -1 for the symmetry elements listed in the caption 
of Fig. 7. The magnetic moments of Ni2 atoms are restricted by symmetry to be parallel 
to 6-axis, thus forming a ferromagnetic sublattice. The y-components of Nil spins are also 



6 




FIG. 3: (Color online) (a) - (f) Temperature dependence of the magnetization of Li2Ni2Mo30i2 
at low T in various magnetic fields. Open and closed circles indicate data measured in the RFC 
and FC processes, respectively, (g) Temperature dependence of the magnetic susceptibility and its 
inverse in the paramagnetic state. 
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FIG. 4: (Color online) Magnetic field dependence of the magnetization of Li2Ni2Mo30i2 at 1.7 K. 
Arrows indicate directions of magnetic field scans. 
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FIG. 5: (Color online) A neutron powder diffraction pattern of ^Li2Ni2Mo30i2 at 10 K (higher 
than Tc) measured using the HRPT diffractometer (A = 1.886 A). Lines on the observed pattern 
and at the bottom show a Rietveld refined pattern and difference between the observed and the 
Rietveld refined patterns. Hash marks represent positions of nuclear reflections. 

forced to be ferromagnetically aligned, whereas the spins in a;2;-plane have antiferromagnetic 
configuration with basis vectors listed in the caption of Fig. 7. 

In the refinement of the magnetic structure, all the atom structure parameters were fixed 
by the values obtained from HRPT data (Table I). The magnetic structure is depicted in 
Fig. 7. The magnetic moment at the Nil(l) position is (0.27(4), 0.04(12), -1.58(3))/iB/Ni 
and its magnitude is 1.60(3) /xb/NI at 2.0 K . The ordered Nil moment mainly points along 
the c direction with a small component along the a direction. The size of the ferromagnetic 
h component is refined to a small value close to zero within experimental errorbars. The 
arrangement of the moments indicates that the Ji, J2, and J3 interactions are FM, FM, and 
AFM, respectively. The magnetic moment at the Ni2(l) position is (0, -0.6(1), 0)/iB/Ni. 
Ni2 moments form the FM LRO. The Ni2 moment is reduced in comparison with ordinary 
ordered Ni^"*" moments. The reduction may be caused by the one dimensionality of the Ni2 
spin system and large disorder due to about 50 % substitution of Li ions. 

The IR Ts well fits experimental data for all the measured temperatures. Figure 8 shows 
the T dependence of the magnetic moments. The absolute value of the FM component m2b 
increases monotonically with decreasing T. The AFM components mia and mic have signif- 
icantly weaker dependence of temperature with step like change at Tc. The FM component 
mife stays close to zero with slight tendency to increase. 
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TABLE I: Structural parameters of Li2Ni2Mo30i2 derived from Rietveld refinement of the HRPT 
neutron powder diffraction pattern at 10 K. The space group is orthorhombic Pnma (No. 62). The 
lattice constants at 10 K are a = 5.0639(5) A, 6 = 10.397(1) A, and c = 17.467(2) A. Estimated 
standard deviations are shown in parentheses. The atomic displacement parameters B[so were 
constrained to be the same for the same atom types and for the atoms occupying the same site. 
The reliability factors of the refinement amounted to i?wp = 3.21, i?exp = 1.77, = 3.31, and 
-^Bragg = 2.76 
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To further elucidate the static and dynamic magnetic properties of Li2Ni2Mo30i2, we 
performed zero field /iSR experiments for various temperatures. In Fig. [9] zero field /iSR 
spectra are shown for characteristic temperature above and below Tc = 8.0 K. Figures |9] (a) 
and (b) show the same data but plotted on different time scales to highlight the observed 
changes for the fast and slow relaxing sig nals. At high temperatures a slowly relaxing 
spectrum of Gauss-Kubo-Toyabe form 43| is observed. This is the fingerprint of a non- 
magnetic material. The small Gaussian relaxation anm is due to static random nuclear 
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FIG. 6: (Color online) (a) A fragment of the difference between two neutron powder diffraction 
patterns of ^Li2Ni2Mo30i2 at 2.0 K and 10 K (magnetic neutron diffraction pattern) measured 
using the DMC diffractometer (A = 4.206 A) Lines on the observed pattern and at the bottom show 
a Rietveld refined pattern in le Bail profile matching mode and difference between the observed 
and the Rietveld refined patterns. Hash marks represent positions of magnetic reflections. Several 
indices of the reflections are written, (b) A neutron powder-diffraction pattern at 2.0 K. Lines on 
the observed pattern and at the bottom show a Rietveld refined pattern for the magnetic structure 
model Ts and difference between the observed and the Rietveld refined patterns. Upper and lower 
hash marks represent positions of nuclear and magnetic reflections, respectively. 



moments only. Below T^. = 8.0 K an over-damped spontaneous muon spin precession is 
observed, indicating a very broad field distribution at the muon site(s), i.e. a complicated 
or disordered magnetic structure or many muon sites in the crystallographic lattice. This 
result is consistent with structural disorder caused by Li ions. 

The /iSR spectra have been fitted using the following functional form for the polarization 
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FIG. 7: (Color online) The magnetic structure of Li2Ni2Mo30i2. Ni2 moments form FM LRO. 
To show clearly the arrangement of Nil moments, symmetry operators of Ml (8d) sites are given. 
(1) X, y, z[u, V, w], (2) X + 1/2, y,z + l/2[u, v, w], (3) x,y + 1/2, z[n, v, w], (4) x + 1/2, y + 1/2, z + 
l/2[u,v,w], (5) x,y,z[u,v,w], (6) x + 1/2, y, z + l/2[u, tD], (7) x,y + 1/2, z[u,v,w], and (8) 
X + 1/2, y + 1/2, z + l/2[u, v, w]. Here u, v, and w indicate components of the magnetic moment. 



function: 



with 



'^(^) (-^ fmag^ ^nmi,^^ ~l~ fmag ^mag(^) (l) 



and 



G^ag{t) = Q cos(27r^t) e-^^-^--*)' + 1^ e-v^. (3) 

Here Gmag and Gnm represent the relaxation functions for the magnetic {fmag) and non- 
magnetic (1 — fmag) voIumc fractious with their corresponding relaxation rates amag and 
(T„m, respectively. Xdyn is the dynamic relaxation rate. 

Figure [10] shows the obtained parameter values from the fits using the above equations. 
A sharp magnetic transition at = 8.0 K is observed in all the parameters. The magnetic 
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FIG. 8: (Color online) Temperature dependence of the components of ordered magnetic moments 
in Li2Ni2Mo30i2. (a) —rriic and mia- (b) rriih and m2b- (c) The average spontaneous magnetization 
per Ni mg. The solid line represents the RFC magnetization shown in Fig. 3(a). 



fraction f^ag reaches ~ 100 % shortly below T^. The magnetic order parameter measured via 
the /iSR frequency as well as the static field width amag continuously increase below Tc. 
The values of <Jmag are larger than the values of indicating the very broad width of the 
internal field distribution. Approaching the transition from above an additional exponential 
damping which increases towards Tc is observed. This is typical for the slowing down of 
electronic moments approaching Tc. The corresponding dynamical relaxation rate Xdyn is 
shown in Fig. [10] (d). Below the transition A^y^ is nearly zero indicating a static magnetic 
state. Decreasing the temperature further, a dynamic relaxation of the 1/3 tail is observed 
again. This is already clearly visibly from the comparison of the /zSR spectra at 7.6 K 
and 1.9 K in Fig. [9] (b). The dynamic relaxation at low temperature had to be fitted with 
a root-exponential time dependence (Eq. |3]), which is indicative of a broad distribution of 
fluctuation times. Interestingly, the dynamic relaxation rate starts to increase below ^ 4K 
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FIG. 9: Zero field muon spin relaxation spectra for various characteristic temperatures shown on 
a short (a) or a long (b) time scale. 



where the magnetization becomes negative. 



IV. DISCUSSION 



Let us consider the origin of the negative magnetization in Li2Ni2Mo30i2. The Neel model 
requires two FM sublattices that in our case are formed by the mif, and m2b components. 
The values of are negative. The values of rriih are very small and show both the signs. 
At 2.0 K, the error bar of rriib is the smallest and the value of mib is positive. Therefore, 
if the T dependence of rriib is similar to that of —rriic or rriia, mif, may be intrinsically 
positive. Compensation of rriib and m2b is possible. The sublattice magnetizations Ma 
and Mb defined in the Introduction correspond to rriib and m26, respectively. Ma has a 
weaker T dependence than Mb in the Neel model. If rriib is proportional to —niic or rriia, 
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FIG. 10: Parameter values obtained from the fits to tlie zero field muon spin relaxation measure- 
ments. The lines are guides to the eye. 
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mib has a weaker T dependence than m2b as the Neel model. We calculated the average 
spontaneous magnetization per Ni rris = pimib + {l —pi)m2b and plotted it together with the 
RFC magnetization in Fig. 8(c). Here, the proportion coefficient pi = 0.728 is calculated 
from the occupancies listed in Table I. The contribution of rriib is large in mg because of the 
large number of Nil atoms. Although the component m2b shows well defined decrease with 
temperature lowering [Fig. 8(b)], nis is practically constant [Fig. 8(c)]. The difference in 
the T dependence of m2b and rris suggests compensation of the FM components. A reversal 
of rris, however, was not obtained. The absolute values of rus seem much larger than the 
absolute values of the RFC magnetization represented by the solid line in Fig. 8(c). 

The negative magnetization is too small to be detected using neutron powder diffraction 
techniques. Therefore, we could not verify the validity of the Neel model. The magnetic 
scattering length was calculated as 2.7gSf{Q) fm using a Born approximation. Here, g and 
S are the g value and magnitude of ordered spins, respectively. f{Q) is the magnetic form 
factor and its magnitude is in the order of unity or less. The nuclear coherent scattering 
length of Ni was obtained as 10.3 fm (natural average). 4^ The intensity of Bragg peaks is 
proportional to the square of scattering lengths. The lower limit of magnetic moments that 
can be evaluated in neutron powder diffraction experiments is about 0.1 fiB- As is seen in 
Fig. 8, error bars are in the order of 0.1 fiB- 

The following four alternative models may be considered as the candidates for explanation 
of the negative magnetization effect. In model 1 applicable to GdCrOs or (Lai?£')Cr03 {RE 



Pr or Nd), 19l-l24| magnetic moments of Cr ions form a canted AFM LRO, while RE 



ions are paramagnetic. A Dzyaloshinskii-Moriya (DM) interaction exists between Cr and 
RE spins. The main (x) component of the Cr moment produces an internal magnetic field 
parallel to the z direction on RE ions. The polarization of paramagnetic moments on RE ions 
induced by the internal magnetic field is aligned opposite to the canted components parallel 
to the z direction. Magnitude of sum of the canted components is larger than magnitude 
of sum of the polarization of the paramagnetic magnetization below Tn, while the latter 
becomes larger than the former below Tcomp- The direction of the main component of the 
Cr moments is fixed. Therefore, the direction of the canted component is also fixed. As a 
result, the negative magnetization appears in GdCrOa or (Lai?i?)Cr03 {RE = Pr or Nd). 
In Li2Ni2Mo30i2, the distance between Nil and Ni2 sites is larger than 5.19 A. Probably, 
exchange interactions between Nil and Ni2 spins are weak. Internal magnetic fields produced 
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by Nil (Ni2) moments on Ni2 (Nil) sites are small. Accordingly, the negative magnetization 
in Li2Ni2Mo30i2 cannot be explained by the model 1. 

In model 2 applicable to LaVOg 26|-28| or YV03,[29|-l3l| only the canted component of the 



magnetic moment on V ions reverses and becomes anti-parallel to external fields. Therefore, 
the negative magnetization was observed. In LaVOa, a structural phase transition occurs 
just below Tn due to the Jahn- Teller effect. As a result, the DM vector reverses and therefore 
the canted component of the magnetic moment reverses. In YVO3, the canting direction is 
determined mainly by the single-ion anisotropy at high T and by the DM interaction at low 
T. As a result, the canting direction reverses at low T. In Li2Ni2Mo30i2, LRO of Ni2 spins 
is the FM order (not canted AFM order). Consequently, the model 2 cannot account for the 

negative magnetization in Li2Ni2Mo30i2: 

In model 3 applicable to (SmGd)Al2, 32| compensation between the spin and the orbital 
parts of the ordered moments occurs. We observe a sum of spin and orbital magnetic 
moments. Therefore, the results in Fig. 8 indicate directly that compensation between 
the spin and the orbital parts of the ordered moments is not the origin of the negative 
magnetization. In addition, in 3d shells, only spin magnetic moments exist because of 
quenching of the orbital angular momentum. 



In model 4 possibly applicable to Sr2i?Mlu06 {RE= Yb, Y or Lu) . |33l-l36| exchange bias 
(EB) effects [45] may cause negative magnetization. The effects can appear in composites 
having both FM and AFM LROs. Ordered magnetic moments of the two LROs should 
be collinear. The FM transition temperature should be higher than the AFM transition 
temperature. In Li2Ni2Mo30i2, however. Nil and Ni2 moments are nearly perpendicular to 
each other. Only one transition was observed within experimental accuracy. In addition, 
distances in Nil - Ni2 bonds are larger than 5.19 A, suggesting that interactions between 
Nil and Ni2 spins are weak. Consequently, EB effects cannot account for the negative 
magnetization in Li2Ni2Mo30i2. 

At present, we think that the Neel model may be valid for the negative magnetization in 
Li2Ni2Mo30i2, although we cannot prove the validity of the Neel model. We cannot deny 
perfectly possibility of another (unknown) origin. As was described in the yuSR results, the 
dynamic relaxation rate starts to increase below ^ 4 K where the magnetization becomes 
negative. It is tempting to conclude that the magnetization inversion mechanism involves a 
dynamical process. As another idea, a structural change, which is not detected, may cause 
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reversal of m2b- Experimental investigations into Li2Ni2Mo30i2 using other techniques and 
further theoretical considerations are necessary. 

V. SUMMARY 

Wc studied magnetism of the spin-1 substance Li2Ni2Mo30i2. The spin system is com- 
posed of the distorted honeycomb lattices and linear chains formed by two distinct Ni^"*" 
sites, Nil and Ni2 sites, respectively. Li+ ions substitute about 25 % and 50 % of the hon- 
eycomb and chain Ni sites, respectively, creating the disorder in both the spin subsystems. 
A magnetic phase transition occurs at Tc — 8.0 K. In low magnetic fields, the magnetiza- 
tion increases rapidly just below Tc, decreases below 7 K, and finally becomes negative at 
low temperatures on cooling. We determined the magnetic structure using neutron powder 
diffraction results. The magnetic order on both Ni-sites develops according to a single ir- 
reducible representation T5 of the space group Pnma and propagation vector k = 0. The 
honeycomb lattices show antiferromagnetic (AFM) long-range order. The ordered Nil mo- 
ment mainly points along the c direction with a small AFM component along the a direction. 
The size of the ferromagnetic (FM) b component is very small. The ordered Ni2 moment 
on linear chains has only the b component and shows FM long-range order. The magnetic 
moment sizes amounted to 1.60(3)and 0.6(1) hb pcr Ni-atom for the honeycomb-lattice and 
chain Ni sites, respectively. The sharp magnetic transition at Tc was also observed in the 
//SR results. The over-damped spontaneous muon spin precession below Tc indicates a very 
broad field distribution at the muon site(s). The dynamic relaxation rate starts to increase 
below 4 K where the magnetization becomes negative. We have examined several models 
known from literature that can account for the negative magnetization. The Neel model 
that is based on the compensation of the ferromagnetic sublattices may be valid in our case, 
but our experimental data does not allow us to prove this model quantitatively. 
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